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Abstract 


We compare the capability of wavelet functions used for noise removal in preprocessing step 
of a QRS detection algorithm in the electrocardiogram (ECG) signal. The QRS signal to noise 
ratio enhancement and the detection accuracy of each wavelet function are evaluated using three 
measures: (1) the ratio of the maximum beat amplitude to the minimum beat amplitude (RMM), 
(2) the mean of absolute of time error (MATE), and (3) the figure of merit (FOM). Three wavelet 
functions from previous well-known publications are explored, i.e., Biorl.3, DblO, and Mexican 
hat wavelet functions. Results evaluated with the ECG signal from MIT-BIH arrhythmia database 
show that the Mexican hat wavelet function is better than the others. While the scale 8 of Mexican 
hat wavelet function can provide the best enhancement in QRS signal to noise ratio, the scale 4 
of Mexican hat wavelet function can provide the best detection accuracy. These results may be 
combined and may enable the use of a single fixed threshold for all ECG records leading to the 
reduction in computational complexity of the QRS detection algorithm. 
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I. INTRODUCTION 


The detection of the heart disease is very important. One of the measurements used for 
checking the condition of the heart is the electrocardiogram (EGG) signal. The components 
of each EGG beat consist of the P wave, the QSR complex, and the S wave. However, 
the manual analysis of EGG signal is tedious and time consuming. Therefore, a computer- 
assisted system for analyzing EGG signals has been developed to solve these problems. The 
detection algorithm for R peak in the QRS complex is a preliminary step in the computer- 
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and heart 


assisted system to other EGG signal analysis such as arrhythmia analysis 
rate variability (HRV) analysis |3|. 

The algorithm for detecting R peak in the QRS complex is generally divided into three 
main steps, i.e. noise removal, envelope signal determination, and R peak detection. In noise 
removal step, a variety of noises in the EGG signal, such as the muscle noise, the power line 
noise, and the baseline drift noise are removed, one of the most popular methods used for 


inique including continuous 
8l4l4|. and wavelet package 


noise removal in the EGG s^nal is the wavelet transform tec; 
wavelet transform (GWT) discrete wavelet transform 

decomposition 1^. After noise removal, the envelope signal is obtained from the EGG signal 
after noise removal and the threshold is defined to determine the time period where the R 
peak in QRS complex locates. Generally, there are four types of thresholding techniques, 
i.e., a single level fixed^hreshold j7[j_^multiple levels fixed threshold 
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adaptive threshold 


121 Il5| , a single level 


8, 9, [ 13 , 14 1, and multiple levels adaptive threshold. The R peak 
can be detected from the time position where the amplitude of the EGG signal after noise 
removal is maximal. 

Generally, more complicated threshold techniques such as multiple levels hxed threshold 
and the single level adaptive threshold are used because the EGG signals after noise removal 
from different beat types, such as normal beat and premature ventricular contraction beat, 
have signihcant difference in their amplitudes due to their difference in frequency compo¬ 
nents. All wavelet functions in previous publications were selected based on their efficiency 
in removing noise. However, their effects on QRS signal to noise ratio and the detection 
accuracy are not yet considered. Therefore, the proposed comparison measurements of the 
wavelet function in GWT are presented in this paper to demonstrate both QRS signal to 
noise ratio enhancement and detection accuracy so that the reduction in the computational 
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complexity of the algorithm in the R peak detection step can be obtained. 

The rest of this paper is organized as follows. Section 2 presents the CWT theory. Section 
3 describes the QRS detection algorithm, the performance measures for the selection of 
wavelet functions, as well as the measures for detection accuracy. Results and discussion are 
given in Section 4. Finally, conclusions are drawn in Section 5. 


II. CONTINUOUS WAVELET TRANSFORM 

The CWT has been gained popular uses for decomposing signals in many applications 
including noise removal in ECG signals. Given the input signal x{t), which is the ECG 
signal in this paper, the CWT of x{t) can be expressed as 


Ta,b — 




t — b 


dt, 


( 1 ) 


where Ta^ is the CWT of x{t), a is the dilation or scale parameter, b is the location 
parameter, and is the complex conjugate of the wavelet function. The scale and the 

wavelet function are two important parameters affecting the performance of noise removal in 
ECG signals. There are various types of the wavelet functions that are successfully used for 
removing noises in ECG signals from previous publications including Biorl.3 0, DblO Q, 


and Mexican hat wavelet 6(] functions. This paper proposes to study and quantitatively 
compare the performance of these three wavelet functions on their capability in enhanc¬ 
ing QRS signal to noise ratio and demonstrate its potential application in QRS detection 
algorithm of ECG data. 


III. MATERIALS AND METHODS 

A. QRS Detection Algorithm 

Fig. [1] shows a block diagram of the QRS detection algorithm used for evaluating wavelet 
functions, which consists of three steps: noise removal using CWT, envelope signal determi¬ 
nation, and R peak detection. Details of each step are as follows. 

1. ) Determine the signal after noise removal y[n] from the ECG signal x[n] by processing 
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FIG. 1: Block diagram of the QRS detection algorithm used for evaluating wavelet functions, 
based on the CWT, which can be expressed as 

yin] = Tl,. (2) 

2. ) Determine the envelope signal z[n] from y[n] using the maximum filter with the 
length L = 120 ms as given by 

z[n] = max y[k]. (3) 

fcE[n—L+l,n] 

3. ) Detect the R peak r[n] in QRS complex using the following steps. 

(a) Dehne a threshold value thv. 

(b) Find the time duration where z[n] is greater than thv and determine the beginning 
time ti and the ending time t 2 - 

(c) Determine the R peak location tn from time between [fi ^ 2 ] in y[n] that gives 
the maximum amplitude. 
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FIG. 2: Signal characteristics in QRS detection algorithm, (a) Thin line: EGG signal after noise 
removal y[n]. Dotted line: Envelope signal 2 ;[n]. Thick line: Threshold value line, (b) EGG signal 
overlaid by the markers from the algorithm (square) and the expert (asterisk). While the signal 
on the left hand side is a normal EGG beat, the signal on the right hand side is a premature 
ventricular contraction beat. 

Fig. [2] shows an example of signal characteristics determined from the QRS detection 
algorithm. Fig. |2](b) shows an example of two beats of the EGG signal x[n\. While the beat 
on the left hand side is a normal EGG beat, the beat on the right hand side is a prematnre 
ventricnlar contraction beat. The noises in the EGG signal x[n\ are removed based on the 
GWT to obtain ?/[n]. Fig. |2](a) shows an example of the EGG signal after noise removal 
y[rt\ as described in the hrst step of the algorithm nsing a thin line. While the normal 
beat has a single peak after noise removal, the prematnre ventricnlar contraction beat has 


4 






















double peaks. Subsequently, the envelope signal z[n\ is calculated from y[n\ as described 
in the second step of the algorithm and is shown in Fig. Wia) using a dotted line. Then, 
the threshold value thv is dehned as shown in Fig. [21^a) using a thick line. Finally, the R 
peak in the QRS complex r[n] is determined from the time where the maximum peak occurs 
within the time period dehned by z[r}\ and the threshold value as described in the third step 
of the algorithm. Fig. Ela) shows an example of the beginning time fi, the ending time ^ 2 , 
and the time where the maximum peak locates. The Zmax and Zmin are maximum beat 
amplitude and minimum beat amplitude, respectively. In addition, the ECG signal overlaid 
by the markers from the proposed algorithm (square) and the expert (asterisk) is shown in 
Fig. [2](b). The time tp and 4 are the position of R peak given by the algorithm and the 
position of R peak given by the expert, respectively. 

B. Performance Evaluation for Wavelet Functions 

Three parameters were used in evaluating the performance of the wavelet function in the 
CWT consisting of the ratio of the maximum beat amplitude Zmax to the minimum beat 
amplitude Zmin (RMM), the sum of absolute of time error (MATE), and the hgure of merit 
(FOM). Their details are as follows. 

• RMM: The RMM can be expressed as 

rMM = (4) 

^min 

If the value of RMM is closed to 1, the amplitude of all beats in the envelope signal z[n] 
is almost equal. In other words, the QRS signal to noise ratio is maximized. As a result, 
the single hxed threshold technique can be used instead of the adaptive threshold 
technique. This is important because it can reduce the computational complexity of 
QRS detection algorithm. 

• MATE: The MATE in the unit of millisecond (ms) is given by 

1 ^ 

MATE = -^|«,(fc)-i,WI. (5) 

^ k=l 

where N is the total number of ECG beats under considerations, tp{k) is the time 
where the R peak locates determined from the algorithm and te{k) is the time where 
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the R peak locates determined from the expert. The best value of MATE is 0, which 
means that all R peaks in the QRS complex are correctly detected by the algorithm. 


• FOM: The FOM can be defined by 


FOM 


1 

RMM X MATE' 


( 6 ) 


FOM is determined using the combination of RMM and MATE. The higher FOM, 
the more appropriate for the wavelet function in removing noises from the EGG signal 
and equalizing the amplitude of EGG beats after processing with the GWT. 


C. Performance Evaluation for Detection Accuracy 


The performance of QRS detection algorithm is evaluated with sensitivity (SEN), positive 
predictive rate (PPR), and detection error rate (DER). SEN is given by 


SEN 


TP 

TP + FN 


X 100, 


(7) 


where true positive (TP) is the number of correct QRS complex predictions. FN is the false 
negative prediction. In other words, the algorithm predicts that there is no QRS complex 
in the location where there is a real QRS complex. PPR can be expressed as 


PPR 


TP 

TP + FP 


X 100, 


( 8 ) 


where FP is the false positive prediction. In other words, the algorithm predicts that there is 
a QRS complex in the location where there is no QRS complex. DER is used for evaluating 
the accuracy of algorithm including both FN and FP values, which can be given by 


DER 


FN + FP 
TP + FN 


X 100. 


(9) 


IV. RESULTS AND DISCUSSION 


A. Comparisons of Wavelet Functions 


We measure the performance of wavelet functions Biorl.3, DblO, and Mexican hat with 
the EGG data record 207 containing the mixture of normal beats and premature ventricular 


contraction beats from MIT-BIH arrhythmia database 


16j. The QRS complex is detected 
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FIG. 3: Comparisons of the RMM value as a function of scale from three wavelet functions. Results 
from the Biorl.3, DblO, and Mexican hat wavelet functions are shown using square, circle, and 
asterisk markers, respectively. 

using the EGG signal from channel 1 or lead II only. In other words, the EGG signal from 
channel 1 was represented by x[n]. 

Fig. inishows the results of RMM as a function of the scale from 1 to 8 applied on the EGG 
signal record 207 from time 12.6 s to 22.6 s. In other words, the GWT processing method 
is when the scale a is varied from 1 to 8. Results from the Biorl.3, DblO, and Mexican 
hat wavelet functions are shown using square, circle, and asterisk markers, respectively. The 
RMM values from all wavelet functions tend to decrease when the scales increase from 1 
to 8. The RMM values from the DblO wavelet function are greater than those from other 
wavelet functions at same scale (Except the scale 1). The maximum RMM is 11.69 at the 
scale 4 of the DblO wavelet function and the minimum RMM is 1.51 at the scale 8 of the 
Mexican hat wavelet function. 

Fig. m shows the results of MATE as a function of the scale from 1 to 8 applied on 
the EGG signal record 207 from time 12.6 s to 22.6 s. Results from the Biorl.3, DblO, 
and Mexican hat wavelet functions are shown using square, circle, and asterisk markers. 
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FIG. 4: Comparisons of the MATE value as a function of scale from three wavelet functions. 
Results from the Biorl.3, DblO, and Mexican hat wavelet functions are shown using square, circle, 
and asterisk markers, respectively. 

respectively. The MATE values from all wavelet functions tend to increase when the scales 
increase from 1 to 8. The MATE values from the DblO wavelet function are greater than 
those from other wavelet functions at same scale (Except the scale 2). The maximum MATE 
is 56.4 ms at the scale 6 of the DblO wavelet function and the minimum MATE is 1.7 ms 
at the scale 4 of the Mexican hat wavelet function. 

Fig. |5] shows the results of FOM as a function of the scale from 1 to 8 applied on the EGG 
signal record 207 from time 12.6 s to 22.6 s. Results from the Biorl.3, DblO, and Mexican 
hat wavelet functions are shown using square, circle, and asterisk markers, respectively. The 
FOM values from hrst Eve scales of the Mexican hat wavelet function are greater than those 
from other wavelet functions at the same scale. However, at the scale 6, 7 and 8, the FOM 
values from the Biorl.3 wavelet function are greater than those from other wavelet functions. 
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FIG. 5: Comparisons of the FOM value as a function of scale from three wavelet functions. Results 
from the Biorl.3, DblO, and Mexican hat wavelet functions are shown using square, circle, and 
asterisk markers, respectively. 

B. Signal Characteristics 

To gain more insight of the performance measurement, signal characteristics from the 
scale 8 of the Mexican hat wavelet function, the scale 4 of the Mexican hat wavelet function, 
and the scale 4 of the DblO wavelet function are shown and discussed. While the signal 
characteristics from the scale 8 and 4 of the Mexican hat wavelet function are the represen¬ 
tation for the best RMM and MATE, respectively, the signal characteristics from the scale 
4 of the DblO wavelet function are the representation of the worst FOM. 

Fig. |6] shows signal characteristics of the proposed algorithm from the scale 8 of the 
Mexican hat wavelet function applied on the EGG signal record 207. Fig. E^a) shows the 
EGG signal before noise removal x[n] from time 12.6 s to 22.6 s. The EGG signal x[n\ 
consists of 5 normal beats and 5 premature ventricular contraction beats. Fig. |6]^b) shows 
the EGG signal after noise removal y[vt\. While the normal beat has a single peak amplitude, 
the premature ventricular contraction beat comprises double peak amplitudes where the left 
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FIG. 6: Signal characteristics from the scale 8 of the Mexican hat wavelet function applied on the 
EGG signal record 207. (a) The EGG signal before noise removal x[n]. (b) The EGG signal after 
noise removal y[n]. (c) The envelope signal z[n]. (d) The EGG signal overlaid by the markers from 
the proposed algorithm (square) and the expert (asterisk). 

peak is greater than the right peak. Fig. [6](c) shows the envelope signal z[n]. While the hfth 
beat has maximum amplitude, the forth beat has minimum amplitude. As a result, the best 
RMM value of 1.51 is determined. Fig. [6](d) shows the EGG signal overlaid by the markers 
from the proposed algorithm (square) and the expert (asterisk). While all R peaks in QRS 
complex of normal beats are correctly detected, all R peaks in QRS complex of premature 
ventricular contraction beats are slightly incorrectly detected because the algorithm detects 
the left peak, which is greater than the right peak in all premature ventricular contraction 
beats. This leads to the MATE value of 34.3 ms. The corresponding FOM value in this case 
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FIG. 7: Signal characteristics from the scale 4 of the Mexican hat wavelet function applied on the 
EGG signal record 207. (a) The EGG signal before noise removal x[n]. (b) The EGG signal after 
noise removal y[n]. (c) The envelope signal z[n]. (d) The EGG signal overlaid by the markers from 
the proposed algorithm (square) and the expert (asterisk). 

is 0.02. 

Fig. [7] shows signal characteristics resulting from the scale 4 of the Mexican hat wavelet 
function applied on the EGG signal record 207. Fig. [71(b) shows the EGG signal after noise 
removal y[n\. Similar to the case in Fig. |6l while the normal beat has single peak amplitude, 
the premature ventricular contraction beat comprises double peak amplitudes. However, the 
position of the larger peak amplitude in each premature ventricular contraction beat is on the 
right hand side. Fig. dc) shows that the maximum beat amplitude and the minimum beat 
amplitude locate at the forth beat and the fifth beat, respectively. Gonsequently, the RMM 
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FIG. 8: Signal characteristics from the scale 4 of the DblO wavelet function applied on the EGG 
signal record 207. (a) The EGG signal before noise removal x[n\. (b) The EGG signal after noise 
removal y[n]. (c) The envelope signal z[n]. (d) The EGG signal overlaid by the markers from the 
proposed algorithm (square) and the expert (asterisk). 

value of 3.28 is obtained. Fig. [71(d) shows the EGG signal overlaid by the markers from the 
proposed algorithm (square) and the expert (asterisk). All R peaks in QRS complex of both 
normal beats and premature ventricular contraction beats are almost correctly detected. All 
R peaks in QRS complex of premature ventricular contraction beats are correctly detected 
because the proposed algorithm detects the right peak, which is greater than the left peak 
in all premature ventricular contraction beats. As a result, the best value of MATE 1.7 ms 
is achieved. In addition, the corresponding value of FOM in this case is the best at 0.18. 

Fig. IHlshows signal characteristics resulting from the scale 4 of the DblO wavelet function 
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applied on the ECG signal record 207. Fig. |8]^b) shows the EGG signal after noise removal 
y[rt\ consisting of multiple peaks, which are different from double peaks in the case of Mexican 
hat wavelet function. Fig. |H]^c) show the envelope signal z[n]. The position of maximum 
z[n\ and the position of minimum z[n\ are at the hfth beat and the eighth beat, respectively. 
This results in the worst RMM value of 11.69. Fig. El^d) shows the EGG signal overlaid by 
the markers from the proposed algorithm (square) and the expert (asterisk). All R peaks 
in QRS complex of both normal beats and premature ventricular contraction beats are 
incorrectly detected. Therefore, the high value of MATE at 41.9 ms is obtained. Moreover, 
the corresponding value of FOM in this case is the worst at 0.002. 


C. Performance Comparisons 


To demonstrate the potential of selected wavelet functions, we apply the QRS detection 


rows the performance comparisons 
8| ), which use the wavelet transform 


algorithm to the whole EGG data of record 207. Ta ble H sf 
of the algorithm with that from other 3 papers ( 1^ [l3| 5 
and single level adaptive thresholding techniques in QRS detection algorithms. Note that 
each R peak from the proposed algorithm is considered as a correct detection when it locates 
within ±50 ms from the position of the R peak given by the expert. The DER value from 
the scale 4 of the Mexican hat wavelet function (1.08%) is smaller than that from the scale 
8 of the Mexican hat wavelet function (1.29%). Results show that the DER values from 
both wavelet functions are in the same range as those from other publications. However, the 
proposed algorithm uses only a single hxed threshold without any additional post-processing 
techniques. 


V. CONCLUSIONS 

We compare the capability of wavelet functions used for noise removal in QRS detection 
algorithm. The effects of the wavelet function on the performance in terms of QRS signal to 
noise ratio enhancement and detection accuracy are carefully studied and evaluated using 
three measurement parameters, i.e., RMM, MATE, and FOM. Three wavelet functions from 
previous publications are explored consisting of the Biorl.3, DblO, and Mexican hat wavelets 
functions. Results show that the Mexican hat wavelet function is the most appropriate for 
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TABLE I: Performance comparisons of the algorithm from scale 4 and 8 of the Mexican hat wavelet 
functions with those from other 3 publications. 


Method 

Total 

TP 

FN 

FP 

SEN(%) 

PPR(%) 

DER(%) 

Scale 4 

1860 

1852 

8 

12 

99.57 

99.36 

1.08 

Scale 8 

1860 

1843 

17 

7 

99.09 

99.62 

1.29 

Choi M 

1860 

1848 

12 

10 

99.35 

99.46 

1.18 

Chen [1^ 

1863 

1860 

3 

24 

99.84 

98.72 

1.45 

Zidelmala [8j 

1872 

1860 

12 

8 

99.36 

99.57 

1.08 


the QRS detection algorithm because it can give good results in terms of both QRS signal 
to noise ratio enhancement (scale 8 of the Mexican hat wavelet function) and detection 
accuracy (scale 4 of the Mexican hat wavelet function), which opens the opportunity for the 
use of a single level fixed threshold for all records of ECG data. 

To achieve a better performance, the results suggest the combination of multiple scale 
wavelet functions instead of using a single scale wavelet function as published in previous 
literature. In other words, the combination of multiple scale wavelet functions may simul¬ 
taneously allow for both the enhancement QRS signal to noise ratio and the increase in 
detection accuracy. 
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